Simulations Reveal Fast Mode Shocks in Magnetic Reconnection Outflows 
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Magnetic reconnection is commonly perceived to drive flow and particle acceleration in flares of 
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solar, stellar, and astrophysical disk coronae but the relative roles of different acceleration mecha- 
nisms in a given reconnection environment are not well understood. While outflow fast mode shocks 
l/^ ■ have been predicted analytically, we show for the first time via direct numerical simulations that 

such shocks do indeed occur in the outfiows of fast reconnection when an obstacle is present. These 
shocks are are distinct from the slow mode Petschek inflow shocks. If Fermi acceleration of electrons 
operates in the weak fast shocks the associated compression ratios will induce a Fermi acceleration 
' particle spectrum that is signiflcantly steeper than strong fast shocks commonly studied, but consis- 

tent with the demands of solar flares. While this is not the only acceleration mechanism operating 
in a reconnection environment, it is plausibly a ubiquitous one. 

PACS codes: 96.60.Iv, 52.35.Vd, 95.30.Qd, 52.35.Tc, 96.60.qe, 98.54.Cm 
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C/5 ■ I. INTRODUCTION 



Magnetic reconnection is a frequently studied process in plasma astrophysics whereby magnetic energy is converted 
into some combination of thermal, non-thermal, and flow kinetic energy. Reconnection underlies models of solar and 
' stellar flares and coronae of astrophysical accretion disks. The process has also been studied in laboratory plasma 
. experiments 

I ' Most work has focused on understanding the reconnection rate and under what circumstances fast reconnection 
occurs. Several distinct paradigms for the onset of fast reconnection in nature have been studied and may operate in 
complementary environments. The first occurs when the ion inertial length is larger than the magnetohydrodynamic 
(MHD) Sweet-Parker resistive length, facilitating current driven plasma instabilities that enhance dissipation near 
the X-point 's'-'s]. A second occurs in the context of turbulent MHD, where the global reconnection rate is enhanced 
^-H ' by the contemporaneous action of many small scale sites 6-8]. A third path to fast reconnection occurs for a long 
^ . enough laminar MHD current sheet that becomes tearing mode unstable and again multiple small scale reconnection 
• ^ ' sites acting together enhance the rate. [oj. 

, There has been less synthesized progress understanding how the converted magnetic energy partitions and on the 
H ■ accelerated particle spectra although different mechanisms have been investigated [l0l - [l6j . Because of the variety of 
' plasma reconnection conditions, it is best to think of reconnection as an acceleration environment rather than a 
single mechanism. Different processes can operate with different strengths depending on the conditions, and whether 
reconnection is fast or slow. 

The basic question we address here, using numerical simulations, is whether fast shocks arise generically in recon- 
nection outflows. Analytic predictions and phenomenological scenarios suggest that weakly compressive fast shocks 
may form downstream of reconnection sites till [l5l - fl7| , such as in coronal flares if supersonic flows from the X-point 
impinge on a previously reconnected loop top from above. However, there has yet to be a study of this specific such 
shock formation numerically. Note that the fast outflow shocks we are focusing on are distinct from slow mode inflow 
shocks of the Petschek model [l^ . Even weak fast shocks in reconnection outflows would be important because of their 
potential ubiquity and role for shock- Fermi acceleration [19, 20]. A fast downstream shock in solar flares [TTI. [TtI. 
also benefits from magnetic dissipation near the X-point that may provide the needed injection electrons that seed 
the Fermi process. 

Here we present 2-D magnetohydrodynamic (MHD) simulations of initial Harris current sheet configurations that 
are unstable to magnetic reconnection, and determine the conditions for which the outflows are weakly super fast- 
magnetosonic. By injecting a dense plasma obstacle into the downstream, we then directly show the formation of 
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weak outflow shocks. From the numericahy computed compression ratios, we discuss the imphcations for electron 
energy spectra from shock Fermi acceleration and discuss how the results may apply to solar flares. 



II. METHODS 



For our simulation we use the MHD code ATHENA [23, l2J|, a cartesian, time explicit, unsplit, Godunov, code 
parallelized by MPI for compressible MHD. The mass conservation, momentum, and magnetic induction equation to 
be used are given by 

|^ + V.(pv) = (1) 
+ V- (pvv - BB) + VP* =0 (2) 

at 

— + Vx(vxB + z/„,J) = 0, (3) 



where P* = P + (B • B)/2, P is fluid pressure, B is the magnetic fleld, p is the mass density, v is the velocity field, 
and the magnetic diffusivity i/,„ equals the resistivity, i.e. Vm = Vm- In the units used in ATHENA the magnetic 
permeability is unity and one converts to cgs units by dividing B by \/47r. 
In Eq. ([3]) J = V X B, and we write Vm as 

V„i = (1^0 + i^ee^), (4) 

where Vc is a uniform background resistivity and i>e is an enhanced resistivity. The quantity a determines the radial 
scale of the enhanced resistivity about the simulation origin surrounding the X-point. The functional form of the 
gaussian fall off of this enhancement is arbitrary and is simply chosen to provide the option to consider the effect of 
enhanced dissipation near the X-point for v,, > 0, the importance of which we will discuss later. 

We set up initial isothermal equilibria for all cases but then distinguish cases in which perturbations are evolved 
adiabatically or isothermally. These cases are referred to as "adiabatic" and "isothermal" respectively. Note that 
by adiabatic we mean no heat loss from the system. Neither the adiabatic or isothermal cases are isentropic. As 
the isothermal and adiabatic simulations bracket a wide range of cooling efficiencies for the perturbed configurations 
motivates our inclusion of both cases. 

For both adiabatic and isothermal simulations we initialize P as 

P = Pcl (5) 

where Cs is the sound speed and is set equal to .5 initially and p is initialized as discussed later in this section. 

For our isothermal simulations no energy equation is solved and Eq. ([5]) becomes a fixed definition. For our 
adiabatic simulations we replace Eq. ([S]) with the equation for energy evolution given by 

OT T 

— = V-[(U + P*)v - B(B • v) + i^^J X B] = 0, (6) 

where 



U = P/(7-l) + (pv.v/2) + (B.B)/2. (7) 

and Cs evolves as yJ^Pj p with 7 = 5/3 and evolves as P evolves according to Eq. (j6l). 

Our initial conditions follow those in the subset of the MHD GEM challenge simulations f?] with only resistivity 
and no viscosity. Technically, in MHD the gyro-radii and plasma particle inertial lengths are infinitesimal, and the 
non-ideal nature of the plasma is scalable entirely with the Lundquist or Magnetic Reynolds numbers given by Lva/vM 
where Va = By^Airp is the Alfven speed. 

It is tempting to normalize lengths to the ion inertial length c/ujpi ~ c/ ATrriQe^ /rm and our time scale to the 
inverse ion gyrofrequency {eB/mic)~^ to see what kind of scales our choices of number density and magnetic field, 
and temperature would imply for a real system as well as to compare to the MHD and particle in cell simulations of 
the GEM Challenge^. In our simulations, x is the inflow axis and y is the axis of the current sheet. We employ a 
rectangular domain with —L^jl < x < and —Ly/2 < y < Ly/2, where L^, = 12.8 and Ly = 25.6 in arbitrary 

units. Upstream conditions of a particular solar flare [l^], no = lO^cm"'^, and Bq — 20G, corresponds to an ion 
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inertial length of 720cm, ion gyroradius of 29cm, and Alfven speed va ^ 1378km s"-'^ so this indicates that that for 
this flare we (and the authors of other simulations would be assuming MHD is valid for scales > 12.8 x 720cm. 
This is misleading however, since the presence or absence of shocks once in the MHD regime does not depend on the 
actual normalization of the density or the magnetic field. One need only be concerned that the MHD approximation 
is valid on scales for which a comparison to a real system is being made and to check the correpsonding Lundquist 
number. 

Resolutions for the runs in Table 1 range from 384 x 768 to 960 x 1920. Most simulations were run at 480 x 960. 
The primary effect of resolution was to change the steepness of the flow profiles near the x point but not the final 
magnitude of the outflow at the simulation edges. Outflow boundary conditions were used at all edges for both fluid 
and magnetic quantities. The mass and magnetic field were not resupplied as the simulations evolved. Perturbation 
around a stable configuration causes the commencement of reconnection, which draws material toward the y-axis. 

ATHENA maintains a zero divergence of the magnetic field if the initial divergence is zero. To ensure V • B = 
in ATHENA, we initialize the magnetic field via a vector potential A such that V x A = B. For our unperturbed 
Harris configuration we used 

= — — log (cosh (a;/a))z, (8) 
a 

where &o = 1-0 and a — 0.5 for all simulations which then returned 

dA 

By{x) — — = 5otanh (a;/a). (9) 

The density profile for our initial Harris configuration is 

Pi^)^ , f ° , , + Poo, (10) 

cosh [x I a) 

where po = 1 and poo = 0.2 for all simulations. Inclusion of a constant background density was used to prevent a large 
an Alfven speed in low density regions that would result in an unreasonably small time step. The overall pressure 
balance and dynamics were left unchanged by the background density. 

Using Cs = a/cTS, makes P* = 6q/2 + PooC^, and /S — 5q/ {2pooC^) = 0.2 and represents an initial equilibrium around 
which we peturb. For the cases we call "adiabatic" we set the initial unperturbed P in Eq. (jH]) to the isothermal 
equilibrium pressure, but then set 7 ~ 5/3. The deviation from equilibrium then evolves adiabatically according to 
Eq. (HD. 

To initiate reconnection, we use a vector potential perturbation 

A'^{x, y) — —ipo cos{kx) cos{ky) (11) 
where k — = and i/jq = 1 for all simulations. The perturbed magnetic field is then 

B' = [sin (kx) cos {ky)Si — cos (kx) sin (fcy)y] . (12) 

Table |T] summarizes selected runs. The A, I, C, and E prefixes indicate adiabatic, isothermal, uniform resistivity 
(ve — 0), and enhanced resistivity (ve ^ 0) respectively. Simulations IS and IL are equivalent to IE, but with half 
and twice the domain size respectively. 

The prefix W indicates simulations in which a dense plasma "wall" of zero temperature is injected at one end of 
the domain perpendicular to the y direction between y — 10.8 and y — 12.8 to facilitate a shock when the outflow is 
super-magnetosonic. The wall is intended simply as a generic obstacle. In a solar flare a natural obstacle arises in 
standard geometries as the downward directed component of the outflow from the X-point impinges upon the closed 



soft X-ray magnetic loop of previously reconnected field lines [24 1. 

The wall is injected only once per simulation at time , but we insert the wall at different times for different runs 
to test the sensitivity to tw The sound speed of wall density is set to zero to prevent the wall from expanding due 
to pressure gradients. The injected wall must supply a density contrast at least as high as that across any shock it 
facilitates. The wall density for the isothermal simulations was twice that of the adiabatic simulations, consistent with 
the fact that the mass conservation jump condition demands a higher compression ratio across isothermal shocks. 
The five wall simulations with enhanced resistivity (y^ > 0) led to stronger outflows than the two wall simulations 
with Vf, —{s. 
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III. RESULTS 



Runs IC, IE, AC, & AE were unstable to reconnection and resulted in accelerated outflows along the x line. 
For these wall-free runs. Table 1 shows that at the final time tj = 75, the uniform {v^ = 0) vs. enhanced (y^jvf. = 10) 
resistivity cases produce Mach numbers (given by Vyj^Jcl + ) of 0.5 and 1.1 for the adiabatic simulations, and 1.1 
and 2.9 for the isothermal simulations. Figure ([1]) and Figure ([2]) show the resulting flow profiles along the x — i) axis 
at times 42 and 50 and demonstrate that the flows had reached a steady state by i = 50. No flows were found to 
occur for test runs in which the perturbation to the magnetic field was not included. 

The adiabatic simulations require I'e > to produce super-magnetosonic outflows but both adiabatic and isothermal 
cases produce higher outflow speeds with > 0. These trends are expected: First, isothermal flows imply a release of 
thermal energy from the system so that the conversion of magnetic to bulk ram pressure in the momentum equation 
is more efficient for isothermal flows, leading to a higher Mach number. Heat is also released at isothermal shocks 
implying higher compression ratios compared to the adiabatic cases. Second, when the central resistivity is enhanced, 
Petschek type "fast" reconnection ensues [1, 01 albeit for the artificial incompressible limit, a centrally enhanced 
resistivity facilitates the rapid rotation of field lines through the plasma near the X-point required for the Petschek 
configuration^. The infiow magnetic field exerts a tension force parallel to the outflow that helps accelerate the 
latter. Slow shocks (not to be confused with the outflow fast shocks) appear in the inflow near the reconnection region 
for all our simulation cases where enhanced resistivity is used. This is shown in Figure ([6]) and Figure (O. Finally, 
lower mach number flows are expected in adiabatic simulations as the temperature of the plasma is increased which 
corresponds to an increase in the sound speed. 

Although our present simulations are strictly MHD, since the ion-inertial length depends only on the plasma density, 
an enhanced resistivity would be expected for fully kinetic versions of our simulated parameter regime based on a 
comparison between the Sweet-Parker thickness and the ion inertial length using such parameters: For example, 
at i = 50 of a typical (e.g. isothermal) simulation we calculate the Sweet-Parker thickness {5sp = L/S^/"^) where 
S = Lva/vm Using L = 12.8, inflow va ~ 1.5 (and /3 = c^/v^ ~ 0.2) and Vc = 0.01 we flnd Sgp ^ 0.29 < = 1 in 
our units. Plasmas with Sgp < (such as pre-flare solar loops [1^), are unstable to enhanced dissipation near the 
X-point i-lill. 

Fig. |3] shows the vertical flow proflle at i = 50 for an isothermal EOS with varying combinations of v^, and 
without the wall. Isothermal simulations were used as the lack of heating allowed for a clearer study of the effects 
of resistivity on the flow profile. The peak outfiow speed toward the edge of the box and the compression ratios are 
quite insensitive to for Ve > 2.5i^c, and insensitive to Vc for Ve = 0. 

Our main result is the formation of shocks in the wall simulations (those with W in Table [l|, across which we 
measured compression ratios r a,t t = 40,45, 50 along the y-axis. Fig. |4] shows snapshots of the vertical flow along 
a; = at t = 50 for adiabatic > cases in which the wall was inserted at times, tw = 0,30,43.5 (cases AEWO, 
AEW30, & AW43 in Table H]) along with an adiabatic i^e = case with wall inserted at t = 30 (ACW30). Only 
in cases with Vg^jvc > 2.5, did a shock form, and in all such cases maintained a compression ratio of ~ 2. The final 
strength of the shock is relatively insensitive to the time the wall was inserted and the shock location converged by 
t = 50 in all cases. 

Fig. [S]is analogous to Fig. 2] for the analgous isothermal runs. For i^e > with walls inserted at = and 30, a 
shock formed with compression ratio ~ 5.8 < r < 6.8 at t = 45 and ~ 5.6 at t = 50. 



IV. POTENTIAL IMPLICATIONS FOR SHOCK-FERMI ACCELERATION AND FLARES 



In the limit that particle velocities are large compared to shock velocities and both are non-relativistic, the ana- 
lytically computed accelerated particle energy spectrum from shock- Fermi acceleration is [l^ ]^{E) oc E^^, where 

2, so s ~ 2.5. (Note that a different proposed first order Fermi acceleration 



_ 1-+1/2 



r-1 



mechanism that also gives s = 2.5 was proposed in jl3| . There, in the context of relativistic particle acceleration in 
active galactic nuclei, Fermi acceleration is argued to occur as the oncoming reconnection inflow streams flow toward 
each other.) 

For a solar flare, the index s at emission sites is related to the photon number spectral index a via s — a + 2, 
(where s — 1 is the electron number spectral index) based on emission from electron-ion interactions at magnetic loop 
foot-points 28 1 . Electron energy spectral indices inferred from typical hard phases of flares [1^ [29'| are 2 ^ s ^ 3.5, 
consistent with the adiabatic prediction of s = 2.5 from our simulations and consistent with scenarios that invoke fast 
shocks in reconnection down-flows ITtI [2H . 

A typical flow travel time between the inferred X-point height and the loop top is ~ 75 seconds [l^l, and the time 
scales for conductive losses [2^] and bremsstrahlung cooling ^ for > 50keV electrons are comparable or within an 
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order of magnitude less. Any cooling that does incur would make s inferred at the emission site an upper limit to s 
at the acceleration site. Note that the shock itself would be expected to be adiabatic because the shock thickness, 
determined by ion gyro-radii, is much thinner than the cooling scale. 

Having chosen initial parameters to match the GEM challenge^ for consistency, we note that the distance bewteen 
the X-point and loop top obstacle in a typical flare is > 10^ times larger than our simulation box given the box density 
if we were to artificially scale dimensions into units of the ion-inertial length based on our chosen density. However, 
our results are strictly MHD and so the fact that the real solar flare is much larger in such dimensionless units makes 
our conclusions that shocks are expected even more robustly relevant for the solar flare scales since the existence of 
shocks does not depend on the density as an independent parameter once in the MHD regime. 

In addition to the mere presence of shocks and the expected power law index, it is noteworthy that the fast recon- 
nection shocks in our simulations extend across the entire outflow width with very little deviation in the compression 
ratio across along the shock (see Fig. 6). All of the inflow then passes through the shock on its way out via mass 
conservation. The efficiency of electron acceleration at the shock would then not be a question of geometric limitation 
but a question of how efficiently the specific type of fast shocks accelerate electrons. This contrasts, for example, 
direct electric field electron acceleration along the x-line, which encounters only a limited electron throughput. A 
detailed study of the microphysics of shock acceleration at specifically weak fast shocks is an ongoing PIC-simulation 
research topic that is beyond the scope of the present work. The extent to which electrons must be "pre-injected" 
above a certain energy threshold 20] to be accelerated at such shocks is not entirely clear. More on this point below. 

Observationally, solar fiares have different populations of electrons and different paradigms have been proposed 
for the potential role that down stream fast shock acceleration might play [ll|, [l^, [13, HH . Some fiares have more 
non-thermal electrons by fraction than others [s^, HJ. In compact-impulsive flares, which have a large percentage of 
non-thermal electrons, there is a knee separating two power laws in the energy spectrum [3^ with the two power laws 
seemingly requiring different mechanisms [2l| : For the particles within the higher energy steeper power law, lower 
energy particles lag behind the higher energy particles, whereas, within the lower energy power law electrons, the 
high energy particles lag the lower energy particles (s^, [sBj . In the shock reprocessing paradigm [2l| , the lower energy 
power law is taken to be the result of stochastic fermi acceleration in turbulent outflows from the reconnection X-point 
and and the highest energy particles are drawn from those which subsequently get injected into the shock as the flow 
passes through. The "obstacle" causing the shock is a reconnected loop in the standard inverse-y picture of solar flares 
but the fraction of particles passing through the shock in the outflow region then depends the convexity of the loop; 
any pre-injection acceleration that might be needed for the particles to engage in the shock-Fermi acceleration could 
be accomplished by the stochastic acceleration region. In principle this can accommodate[2l| the needed ^ lO'^^/s 
rate of electron throughput in typical flares . 

As we stated in the introduction, reconnection is best thought of as an acceleration environment and the relative 
importance of various mechanisms being possibly different in different flares. The discussion of the previous paragraph 
is one example of how such shocks may fit into such a paradigm. However, the main point of our paper here is to 
highlight the robust appearance of such shocks in simulations and their generic presence, regardless of whatever else 
may also be operating to produce observed phenomenology. 

V. CONCLUSIONS 

From 2-D MHD simulations, we have found that super-magneotsonic outflows are a generic feature of fast magnetic 
reconnection simulations and produce shocks upon encountering a plasma wall. We used the same initial parameters as 
those of the GEM challenge MHD simulations to study this question. We find the fast shocks in our adiabatic cases 
(the most relevant for solar flares) to be weakly compressive with r < 2.5. Our results are generally consistent with 
previous analytic predictions of the existence of weak fast mode shocks in reconnection outflows fill] and consistent 
the use of such shocks for particle acceleration in solar flares [l^ [13, [llj . The potential importance of such fast mode 
shocks and the complexity of magnetic reconnection warranted our explicit demonstration of their formation in a 
reconnection numerical simulation. It is also important to emphasize that the weak fast mode downstream shocks are 
distinct from the Petschek slow mode upstream shocks. 

If Fermi acceleration is operative at such shocks, it will produce steeper particle energy spectra than maximally 
compressive r — 4 adiabatic shocks. The resulting steeper particle spectral indices are not inconsistent with inferred 
distributions of solar flare electrons. An important next step is to directly study the microphysics of particle ac- 
celeration at these weakly compressive fast mode reconnection outflow shocks; most work has focused on particle 
acceleration at strong fast mode shocks. 

We have not considered a guide field along the z-axis in the third dimension and depending on on the spatial 
gradient of of such a field the Mach numbers of the outflow may change. The influence of a guide field would likely 
depend on its initial spatial gradient as to whether its net effect is to enhance or reduce the outflow pressure gradient. 
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For example, if the guide field is very stronger at the X-point and weaker farther out, it could increase the outflow 
Mach number. If the field is uniform it would not produce a strong force. This is another topic for future work. 
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TABLE I: Simulation runs. Column 2 is the resolution with symbols 7 = 240 x 480, II = 384 x 768, III = 480 x 960, 
IV = 960 X 1920. The is the final simulation time, Mp is the peak Mach number at the measured time, r is the compression 
ratio across the shocks at the times indicated, and is relevant only for the Mp > 1 simulations in which a wall is injected at 
time t^. NA indicates not applicable as no wall was present to generate a shock. Magnetic diffusivity parameters i^c, t'e and a 
are defined in the text. 



Sim. # 


res. 






a 




if 


Mp 


r{t = 40/45/50) 


AC 


// 


.01 


0.0 


NA 


NA 


75 


0.46 


NA 


AE 


// 


.01 


10.0 


0.25 


NA 


75 


1.07 


NA 


AEWO 


III 


.008 


10.0 


0.25 


0.0 


50 


1.85 


1.9/2.0/2.2 


AEW30 


HI 


.008 


10.0 


0.25 


30.0 


50 


1.97 


1.8/1.8/2.1 


AEW43.5 


HI 


.008 


10.0 


0.25 


43.5 


50 


1.79 


NA/2.1/2.1 


ACW30 


III 


.008 


0.0 


0.25 


30.0 


50 


0.35 


0.0/0.0/0.0 


IC 


II 


.01 


0.0 


NA 


NA 


75 


1.08 


NA 


IE 


II 


.01 


10.0 


0.25 


NA 


75 


2.91 


NA 


lEWO 


III 


.008 


10.0 


0.25 


0.0 


50 


2.55 


6.0/5.8/5.5 


IEW30 


III 


.008 


10.0 


0.25 


30.0 


50 


2.6 


7.0/6.8/5.7 


IL 


IV 


.008 


10 


0.25 


NA 


50 


2.86 


NA 


IS 


I 


.008 


10 


0.25 


NA 


50 


2.85 


NA 


IDCO 


III 


.016 








NA 


50 


1.1 


NA 


IDC 


III 


.016 


10 


0.25 


NA 


50 


2.92 


NA 


IDE 


III 


.008 


20 


0.5 


NA 


50 


2.92 


NA 


IDEC 


III 


.016 


20 


0.5 


NA 


50 


2.83 


NA 


IHEC 


IV 


.008 


2.5 


.125 


NA 


50 


2.83 


NA 




FIG. 1: Mach number M = Vyj ^ (ci + v\) at a; = for the adiabatic simulations AC and AE at times 42 and 50. Constant 
resistivity simulation AC produced subsonic flows. Simulation AE using enhanced resistivity produced moderately supersonic 
flows. The line marked 'flducial' represents a non-perturbed simulation which remained stable against reconnection and hence 
resulted in no flow velocty. 




FIG. 2: Mach number M = Vyj \J (ci + f^) at a; = for the isothermal simulations IC and IE at times 42 and 50. Constant 
resistivity sinmlation IC produced marginally supersonic flows. Simulation IE using enhanced resistivity resulted ing reatly 
increased flow velocity. One again, the non perturbed simulation marked as 'fiducial' showed stability against reconnection and 
no resulting flow. 




FIG. 3: Mach number M = Vy/^J{ Cs + f a) at X = for the isothermal simulations with varying resistivities. The diagonal line 
represents uniform resistivity i^c > with Vc = 0. Doubling Vc with Ve = leads to nearly the same line. In contrast, including 
enhanced resistivity Ve > Q with strengths ranging from 25% to 200% of Vc, produce the curved profiles which become nearly 
equivalent at large y. 
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FIG. 4: Mach number M = Vy/^/( ci +v\) at X = and t = 50 for the adiabatic simulations of Table 1 with walls placed 
near y 11 at times 0, 30, and 43.5. At negative y there is no wall. The curve that remains subsonic at all y is run ACW30 
in Table 1, the only case shown with = 0. 




FIG. 5: Same as Fig. 2 for the isothermal simulations with walls in Table 1. 
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FIG. 6: Snapshot of the magenetosonic mach number for IEW30 (top image) and AEW30 (bottom image) at t = 50. The 
color bars denote the magnitude of the mach number. The location of the Petschek shocks and the fast shocks are labeled for 
the top figure only but occur at the same location for in the bottom figure, (note the axes in this caption have been flipped to 
allow for easy plotting). 
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FIG. 7: Snapshot of the density profile for IEW30 (top image) and AEW30 (bottom image) at t = 50. The color bars denote 
the magnitude of the density. 



